Static and dynamic properties of vortices in anisotropic magnetic disks 
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We investigate the effect of the magnetic anisotropy {Kz) on the static and dynamic properties 
of magnetic vortices in smaU disks. Our micromagnetic calculations reveal that for a range of Kz 
there is an enlargement of the vortex core. We analyze the influence of Kz on the dynamics of 
the vortex core magnetization reversal under the excitation of a pulsed field. The presence of Kz^ 
which leads to better resolved vortex structures, allows us to discuss in more details the role played 
by the in-plane and perpendicular components of the gyrotropic field during the vortex-antivortex 
nucleation and annihilation. 
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The manipulation of magnetization in nanostructured 
materials by means of magnetic field and/or spin polar- 
ized current has attracted substantial attention in the 
last decade. More recently, a special focus on the study 
of magnetization reversal dynamics in magnetic disks [l[ 
and lines [2] has been motivated mainly by their potential 
importance in the implementation of memory and logical 
operations. 

Micro-sized Permalloy (Py, FeNi) disks can exhibit, de- 
pending on their size and aspect ratio, a magnetic vortex 
with a core (~ 10 — 20 nm) magnetized perpendicular 
to the disk plane [3]. Due to the magnetic bi-stability 
of the vortex structure, they have been considered for 
technological applications. The vortex core magnetiza- 
tion reversal can be achieved by applying an in-plane 
magnetic field or spin polarized current in the form of 
shortpulses [j and/or alternating (AC) resonant excita- 
tion [il, m, - both have an equivalent role on determining 
the vortex core dynamics. 

It is usually considered that the size of the vortex core 
depends on parameters such as exchange constant, thick- 
ness and diameter of the magnetic disk. Most of the 
research on vortices in magnetic systems neglects the ef- 
fect of magnetic anisotropy. On the other hand, it has 
been demonstrated that a uniaxial magnetic anisotropy 
in Permalloy particles can be induced by the deposi- 
tion process [7]. Experiments and simulations have also 
shown that the presence of anisotropy in thick magnetic 
disks gives rise to a diversity of domain patterns ^7|, i8j] . 

In this letter, we study the role played by the mag- 
netic anisotropy in the magnetic properties of disks. Us- 
ing micromagnetic simulations, we consider a magnetic 
anisotropy {Kz) perpendicular to the disk plane and an- 
alyze how it modifies the magnetization pattern and the 
dynamics of the vortex core under action of an in-plane 
pulsed magnetic field. We also discuss in detail the core 
reversal process and the influence of Kz on the gyrotropic 
field and the magnetization reversal time. 

The simulations were performed with a code we devel- 
oped, which employs the Landau-Lifshitz-Gilbert (LLG) 
equation. We used the typical magnetic parameters of 



Py: the saturation magnetization is given by Ms — 
8.6 X 10^ A/m, the exchange coupling is A = 1.3 x 10~^^ 
J/m and the Gilbert damping constant is a = 0.2. The 
magnetic anisotropy was included in the total effective 
magnetic field ^eff and is given by (2Kzl [i^M^^rrizT^ [9|. 
Kz varies from to 10^ J/m^. We have simulated Py 
disks with the diameter of 300 nm and thickness of 12 
nm. The disk is discretized in cells of 3 x 3 x 3 nm^. 

Let us first consider the vortex core structure in static 
equilibrium (in the absence of an external magnetic field). 
The magnetization pattern of the Py disk presents three 
characteristic regimes, as can be seen in Fig. [T](b-f): for 
Kz below 2.5 x 10^ J/m^ the main consequence of in- 
creasing the anisotropy is an increase of the vortex core 
diameter. The size of the vortex core is measured at half 
of the maximum value of rriz and its dependence with 
Kz is shown in Fig. [TJa). For Kz between 2.5 x 10^ and 
4.0 X 10^ J/m^, rriz exhibits concentric regions with +Ms 
and —Ms, still preserving the core in the center of the 
disk. For Kz between 4 x 10^ and 6 x 10^ J/m^ the core 
disappears and the number of concentric rings increases, 
which also happens for larger disks. Moreover, in these 
two ranges, the in-plane magnetization still preserves the 
vorticity (see Figs(T]i-e). The similarity between these 
patterns and previous observed patterns in Co nanomag- 
nets [10] and thick NiFe nanodisks [7] is noteworthy. In 
both cases, the concentric rings are clearly seen in the 
images obtained by Magnetic Force Microscopy (MFM). 
For Kz above 6 x 10^ J/m^, the vortex structure is lost 
and we see what amount to a single domain in rriz (Fig. 
[If). For larger and thinner disks, the magnetization pat- 
tern is composed by stripes. 

The prospect of having large vortex cores has clear ad- 
vantages for magnetization detection. Still, in order to 
consider the practical aspects it is necessary to study 
the stability of these vortex cores and the possibility 
of switching their magnetization. For this purpose, we 
studied the dynamics of the system under short in-plane 
magnetic field pulses. We simulated the magnetic rever- 
sal process and constructed a switching diagram. Simi- 
lar diagrams for pulse parameters have been constructed 
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FIG. 1: (a) The diameter of the vortex core as a function of Kz 
for a disk with diameter of 300 nm and thickness of 12 nm. 
Panels (b-f) show the magnetization pattern for increasing 
values of Kz^ illustrating the different regimes of (a). The 
colors represent the direction of the in-plane component of the 
magnetization vfixy A vortex structure is given by a clockwise 
color sequence of blue- green-red. 



recently for magnetic disks [^. They show that the op- 
erating field range is narrow: low fields do not produce 
core switching. On the other hand, higher fields can give 
rise to multiple switchings. 
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FIG. 2: (a) Switching vortex core magnetization diagram for 
magnetic field pulse strength and Kz. The pulse duration 
is fixed (263 ps). Red, green, and blue colors represent no, 
one, and multiple switches, respectively. In the grey area the 
vortex core is expelled from the disk, (b) The switching time 
as a function of Kz for a fixed field strength Bo — 64 mT. 

With the intention of discussing the influence of K^ 
on this process, we built the diagram sketched in Fig. 
[2^. We use a pulse of Gaussian form with a fixed pulse 
duration (263 ps) and a variable field strength B^. For 
each value of K^^ we count the number of core magne- 
tization inversions during a single pulse length. Fig. [2^ 
shows three different dynamical regimes in response to 
the exciting field. In the absence of magnetic anisotropy, 
the field strength that is necessary to switch the core 
magnetization is B^ = 60 mT. However, for fields higher 
than 95 mT undesirable multiple switches are produced. 
For disks with i^^ 7^ 0, the pulse necessary to induce 
the switching process has a strength comparable to the 
Kz = case. As can be seen in Fig. [2^, an important 



consequence of increasing Kz is the decrease of the mini- 
mum field necessary to produce a single switch. Such de- 
pendence on Kz, for the minimal field for the switching, 
opens the possibility of producing selective vortex inver- 
sions in a group of magnetic disks with different i^^'s. 
Note that under the action of a pulsed field, the core 
moves towards the disk border in a curved (or spiraling) 
trajectory. When the pulse is over, the core returns back 
to the disk center. Depending on the pulse strength and 
duration, the vortex core can be expelled from the disk 
during this process. Such behavior is not shown in our 
diagram and can be avoided by increasing the disk diam- 
eter. 

Another interesting aspect of the influence of Kz can 
be seen in the switching time Tg. Tg is the interval be- 
tween the pulse start and the complete reversal of the 
vortex core magnetization. We have calculated r^ for a 
given pulse strength {Bq = 64 mT) in the single switch- 
ing regime (see the horizontal line in the Fig|2^). Tg for 
Kz = is 213 ps and it increases monotonically with 
Kz, reaching values c^ 1.4 times larger than Ts{Kz = 0) 
(Fig 12)3). For thinner disks (6nm) Tg can reach up to 
2.5r,(i^, =0). 
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FIG. 3: (a) Time evolution of the gyrotropic field (perpen- 
dicular component). Snapshots of vortices and antivortex for 
(b) Kz = at the maxima of the gyrotropic field and (c) 
Kz = 2.5 X 10^ J/m^ and t = 263 ps. The color map repre- 
sents the perpendicular component of the magnetization (niz) 
and the arrows show the in-plane component. The inset of 
panel (a) gives a transverse view of ruz- 



A closer look at the core reversal dynamics for disks 
with Kz y^ shows that the intermediate processes lead- 
ing to the switching are similar to the ones obtained 
in previous analysis for Kz = [l|, 14, [H, [ll|, [l2|: the 
core shape changes during the magnetization reversal. 



through the formation of the adjacent V-AV (vortex- 
antivortex) pair and the subsequent V-AV annihilation 
and nucleation of the reversed vortex. However, in op- 
position to what is normahy observed in micromagnetic 
calculations, here, due to the increase of the vortex core 
diameter, the formation of the adjacent V-AV pair dur- 
ing the reversal process can be well resolved (see inset of 
Fig. [3] and Figs. [3Kb) and (c)). 




FIG. 4: In-plane component of the gyrotropic field (arrows) 
just before (a) and after (b-d) the nucleation and separation 
of the vortex and antivortex with negative magnetizations. 
The color map represents the out-of-plane magnetization. 

The vortex core magnetization reversal, its switching 
time, the core motion and its deformations can be under- 
stood in terms of the gyrotropic field [l|, UjJ , which acts 
on the vortex core only during its movement. The gy- 
rotropic field is given by h^ = -m x [v- Vm] [l|, [ll| , where 
7 and v are the gyromagnetic factor and the core veloc- 
ity, respectively. Indeed, using Thiele's equation [13] in 
its original form, h^ can be simplified using m instead of 
V, and it is reduced to h^ = -m x m. All calculations 
shown below were performed with this last expression of 

Previous discussions on the gyrotropic field have em- 
phasized the role of the perpendicular component of this 
effective field in the reversal process [llj. To further in- 
vestigate this, we have calculated the perpendicular h^ 
and the in-plane h^^ component as well. First, we discuss 
the time evolution of the perpendicular component h^ for 
Kz=0 and 2.5 x lO^J/m^, which is shown Fig.[3fei. Fol- 
lowing the application of the field pulse, the core moves 
and hg acts on the core at the side opposite to the move- 
ment direction leading to the formation of a peak with 
negative magnetization. At this stage, it is just a peak 
without a vortex structure, h^ increases (in modulus) up 
to the point where the peak is so wide that it leads to 
the nucleation of a V-AV pair, labelled as V~-AV~, re- 
spectively. In conjunction with this V~-AV~ nucleation 
there is a decrease (in modulus) of hg. Subsequently, the 
AV~-V+ annihilation takes place at the h^ divergence. 
This divergence is in agreement with Ref. lU,. Such be- 



havior is independent of Kz but it is worth noticing that 
the separation of the V~-AV~ is better resolved spatially 
and takes place after an interval that is longer than the 
one for K^ = 0. This can be seen in the snapshots in the 
inset of Fig. [3l and in the comparison of Figs. [3] (b) and 
(c). On the right one sees the V+ with a close neighbor, 
at its left, which is the AV~. Further on the left, one 
sees the V~. As can be gathered from the figures, h^ de- 
creases for increasing values of Kz and that is the reason 
for the increase of r^ . 

One of the main advantages of considering Kz 7^ in 
our calculations is the large separation between the Vs 
and AVs involved in the switching process. This allows us 
to analyze in more details the structure of the gyrotropic 
field in between these vortex structures. Normally, due to 
the small distance between them, we cannot resolve the 
differences in the gyrotropic field produced by the move- 
ment of vortices and antivortices and the nucleation of 
the V-AV pair cannot be fully understood from this type 
of analysis. For that, we calculated the h^^ component 
during the V~-AV~ nucleation and separation, which are 
fundamental steps in the process of core magnetization 
reversal. Fig. [4^ shows h^^ for an instant just before 
the V~-AV~ nucleation and separation. The color map 
and the arrows illustrate the z magnetization component 
and h^^, respectively. The red circle is the original V+ 
and the blue one is the negative peak. We can see from 
Figs. [4^-b that h^^ is responsible for the transformation 
of this wide negative peak to a V~-AV~ pair and acts 
as a driven force pushing V~. As can be seen Fig. HJd, 
h^^ pushes the AV~ in the direction of V+, producing 
the pair annihilation. In Fig. \^ it is also possible to 
observe the spin waves generated by the AV~-V+ anni- 
hilation. In the last snapshot, we can see the remaining 
V~(see Fig. |4]i). We would like to stress that to our 
knowledge, this is the first analysis that shows explicitly 
the dynamics responsible for both the AV~-AV~ nucle- 
ation and separation and V~-V+ annihilation. For in- 
stance, considering only the perpendicular component of 
h^, together with the vorticity conservation, one is able 
to explain the AV~-V^ annihilation but not the V~-AV~ 
nucleation and separation process. 

In conclusion, we presented a detailed analysis of the 
influence of Kz on static and dynamic properties of mag- 
netic vortices in disks. We showed that increasing values 
of Kz produce a growth of the vortex core. In addition, 
high values of Kz cause a change in the magnetization 
pattern. We then showed, by means of dynamical cal- 
culations using in-plane magnetic field pulses, that both 
in-plane and perpendicular components of the gyrotropic 
field, hg and h^^, contribute and are fundamental to the 
understanding of the vortex core magnetization reversal 
process. 
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